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be possible if neutrinoless double-beta decay will be observed. The required experimental accuracies 
and uncertainty in the calculated nuclear matrix elements of the process are discussed. 
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I. INTRODUCTION 



The discovery of neutrino oscillations in the SuperKamiokande (atmospheric neutrinos) [1], SNO (solar neutrinos) 
0, KamLAND (reactor neutrinos) 3, MINOS Q (accelerator neutrinos) and other neutrino experiments gives us 
compelling evidence that neutrinos possess small masses and flavor neutrino fields are mixed. 

Neutrino flavor states |z/q, > (a = e, /x, r) are connected to the states of neutrinos with masses rrij by the 

following standard mixing relation 

3 

= ^?7aj>i) (a = e, /i,T), (1) 

where U is the 3x3 Pontecorvo-Maki-Nakagawa-Sakata (PMNS) unitary mixing matrix If massive neutrinos are 
Dirac (Majorana) particles in this case the PMNS matrix contains one (three) CP phases. 

The problem of the determination of the CP phases is one of the major challenges in today's neutrino physics. Some 
information about lepton phases could help to solve the problem of the baryon asymmetry of the Universe jGi]. The 
discovery in the Daya Bay , RENO [8] , T2K Q and Double Chooz experiments that the angle ^13 is relatively 
large opened a possibility of the measurement of the Dirac phase i5 in long baseline accelerator experiments. The 
Majorana CP phases can be determined only through the observation of neutrinoless double /3-decay (Oi^/3/3-decay) 

The CP phases enter into the effective Majorana mass [3, [H, defined as 

mi3i3 



Y^Vlm,. (2) 



This quantity depends also on the neutrino oscillation parameters ^13, b^m\^^ Am^^^^j, the lightest neutrino mass 
and the type of the neutrino mass spectrum (normal or inverted) . 

The effective Majorana mass m^^ can be determined in the neutrinoless double beta decay of even-even nuclei 

{A,Z)-^{A,Z + 2) + 2e- (3) 

by converting the 0iyf3/3-dec'dy half-life to |m^^| using calculated nuclear matrix elements (NMEs). 

The goal of this paper is to discuss a possibility for determining the Majorana CP phases from data of 
experiments of the next generation assuming that the 0^'/3 /3-decay will be observed in these experiments. The problem 
of the experimental accuracies and uncertainties of the calculated NMEs will also be addressed. 



II. NEUTRINO OSCILLATIONS AND EFFECTIVE MAJORANA MASS 



It was proved in the experiments with atmospheric, solar, reactor and accelerator neutrinos that flavor neutrinos 
oscillate from one flavor (electron-, muon-, and tau-) to another due to neutrino mixing and non-zero neutrino mass- 
squared differences. All existing neutrino oscillation data (with the exception of the LSND [20|, MiniBooNE 2l|, short 
baseline reactor (2^ and Gallium [2^) anomalies) are perfectly described by the minimal scheme of three- neutrino 
mixing. 

In the case of Dirac neutrinos the unitary 3x3 PMNS neutrino mixing matrix can be parameterized as follows 

(C12C13 C13S12 e~*''si3 \ 

-C23S12 - e*^Ci2Sl3S23 C12C23 - e*''si2Si3S23 C13S23 (4) 
S12S23 - e"^Ci2C23Sl3 -e*''c23Sl2Sl3 " C12S23 C13C23 / 

where Cij = cos^y , = sin^y. 9i2, ^13 and 023 are mixing angles and S is the CP phase. If neutrinos are Majorana 
particles the matrix U in Eq. (|4]) is multiplied by a diagonal phase matrix P = c?m(7(e''-"^/^~'^-', e*'"^/^~*\ 1), which 
contains two additional CP phases ai and a2- 

With the discovery of neutrino oscillations we know: 

• The values of the large mixing angles 612 and 623 ■ The value of the relatively small angle ^13 measured recently 
in the Double Chooz 10], the Daya Bay 7] and RENO [8] reactor neutrino experiments. 
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FIG. 1. (Color online) The CP phase Q2 (left panels) and difference of phases Q21 (right panels) are plotted, respectively, as 
function of the absolute value of the effective majorana mass for the IH and NH of neutrino masses and chosen value of mo. 



• The solar and atmospheric niass-squared difference^ Am^^j^ ~ ^n^i2 ^^'^ ^^atm ~ ^"^23 (Normal spectrum) , 
Am\rj.^,^ = — Am^g (Inverted spectrum). 

We do not know the value of the lightest neutrino mass, the CP phases and the character of the neutrino mass 
spectrum (normal or inverted). 

From the data of the MINOS experiment [4] it was found that Am^^^j = (2.43 ± 0.13)x 10"^ eV^ From the 
analysis of the KamLAND and solar data it was obtained that tan^ 6*12 = 0.452^g'Q|3 From the global fit 
to all data it was inferred that [24] Am^^^^ = {7.65to.lt) x 10"^ eV^ and sin^ 6*23 = O.bOto'^l. Finally from the 
analysis of the Daya Bay[7] and RENO datad one obtains sin226'i3 = 0.092 ± 0.016 (stat) ± 0.005 (syst) and 
sin^26'i3 = 0.103 ±0.013 (stat) ±0.011 (syst), respectively. 

The effective Majorana mass is given by 

\mpp\ = |c?2C?3e*"imi + s^2C?3e'"'"i2 + s?3™3| (5) 

or by the full expression 

|2 44 2,44 2,4 2 
+2c^2'Sl2'^13™l™2 COS (ai — 02) 

+2c22Ci3S]^3mim3 COS ai + 2s]^2Ci3Si3™2TO3 cos a2- (6) 

From this equation it simply follows that the effective Majorana mass depends on the character of the neutrino mass 
spectrum and three unknown parameters: the lightest neutrino mass and two CP phases. 
In the three-neutrino case two mass spectra are currently possible: 

• Normal Spectrum (NS): mi < m2 < m^: Ami2 ^ Am|3. In this case 

7712 = ^ Am2„„ + ml , ms = ^ Am^^^ + Am^^^ + 

with mo — nil. 



^ We use the following definition Am^ • = rrt? — m? 
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Inverted Spectrum (IS), < mi < 1712'- Amfj ^ lA'^ial- We have 



mi = ^Aml^^ + ml, ma = ^ AmJ^^ + Am^^^ + m§ 



with mo = ma 



Here, mo = toi(to3) is the hghtest neutrino mass. 

We will consider the following two neutrino mass hierarchies: 

1. Normal Hierarchy (NH): mi <^ ma ^ m^: 
In this case for the neutrino masses we have 



mi < y/Am^. m2 ~ ^/Am 
Neglecting the negligibly small contribution of mi we find 

\mp,\' - st^c^Am 



2 

SUNT 



ma 



i2 

''ATM' 



COS a2 



2 

SUN 



For the effective Majorana mass we then have the following range of values 



|Sl2Cl3VA^ 



'sun 



'13 



VAr 



l2 

''ATM 



< Imp/sl < Si 



■"SUN 



s?3VA^ 



l2 

''ATM' 



Using the best-fit values of the mass squared differences and the mixing angles we find 

1.5 meV < |m^^| < 3.8 meV 



(7) 



(8) 



2. Inverted Hierarchy (IH): ma ^ mi < ma: 

In the IH scenario ma <^ ^ Am|^„ and mi ~ m2 — Y^Am^^„. We find 



l"^/3/3|'-cla(l-2s?2C?2)Am2 



„„„ ""pp '-lav^ ""i^^i^y"""ATM /n\ 
= 2c^ Am^ 

^''12*12''13^'"'ATM 

where aia = ai — a^- For the absolute value of the effective Majorana mass we find the range 

|cos20i2|c23^Am2^„ < \mfi^\ < c^^^ Aml^^. (10) 

Using the best-fit values of the parameters we find the following range for |m^^| in the case of the IH: 
18 meV < Imppl < 48 meV. 

The absolute value of the neutrino mass can be determined from the precise measurement of the end-point part of 



the /3-spectrum of the tritium [25[ and other /3-decay measurements [26| . Cosmological observations allow to infer the 
sum of the neutrino masses 

3 

m„,„„ = ^mfc. (11) 

k 

The current limits on mooamo depend on the type of observations included in the fit j27|. The CMB primordial gives 
< 1.3 eV, CMB-l-distance < 0.58 eV, galaxy distribution and lensing of galaxies < 0.6 eV. On the other hand the 
largest photometric red shift survey yields < 0.28 eV [2^. It is expected that future cosmological observables will be 
sensitive to m^„j,„„ in the range (0.05-0.1) eV (see, e.g., the recent summary [27]). Thus, the inverted neutrino mass 
hierarchy will be tested by future precision cosmology. In the case of the inverted mass hierarchy we have for the sum 
of the neutrino masses 



m,„,„„ ~ 2^Am\j,^j ~ 100 meV. (12) 

So, there is a chance that future cosmological observation apparently could restrict the value of mo below 10 meV. 

In Fig. [H by exploiting Eq. the Majorana CP phase aa (or difference of phases aai =03 — 0:1) is plotted as 
function of the absolute value of the effective Majorana mass for chosen values of mo and by assuming the NH (IH) of 
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FIG. 2. (Color online) The difference of CP phases 0:21 = 02 — Qi plotted as function of |m/3/3|°'"' for the inverted hierarchy of 
neutrino masses. The current experimental errors of neutrino mixing parameters and mass squared differences are taken into 
account; mo is taken to be in the range of to 10 meV. The blue, red and orange regions correspond to Im^/jl"'" with error 
{app/\mpp\°^'') of 0%, 15% and 25%, respectively. 

neutrino masses. The second phase ai is considered to be arbitrary. We found that for mo within a range of to 10 
meV a2i depends only weakly on |?tt.^^|, if the IH is considered. It is due to the fact that the third term on the right 
hand side of Eq. ([5]) is small in comparison with the first two terms. A different situation occurs in the case of the NH. 
Results depend strongly on mo, where its value cannot possibly be determined from other experiments or cosmology, 
if it is lower than a few meV. Thus, by measuring Oi^/3-decay it will not be possible to obtain a neutrino mass model 
independent information on the value of at least one of the three CP Majorana phases, if there is a normal hierarchy 
of neutrino masses. In some particle physics models there is a possibility to predict all three neutrinos masses. In 
these cases values for one of the CP Majorana phases can be obtained when considering the second phase to be 
arbitrary and by observing the Oi^/3/3-decay. However, the results depend strongly on the value of mo and the type of 
the neutrino mass hierarchy, normal or inverted. It is interesting to note that for mo = 3 meV the minimal value of 
\miip\ does not appear for the case of CP-conservation (see left upper panel of Fig. [T]) but for a2 — 0.797r. 

III. EFFECTIVE MAJORANA MASS AND THE Oi/^/3-DECAY 

Assuming that the Oi^/3/3-decay is driven by the Majorana neutrino mass mechanism, we have for the effective 
Majorana mass mpp 

l"^/3/3| = , „ — • (13) 

Here, Goi^{Qi3i3, Z) and Af'"'^ are, respectively, half-life of the Oi^/3/3-decay, the known phase-space factor and the 
nuclear matrix element, which depends on nuclear structure. 

The future experiments, CUORE (^^^Te), EXO, KamLAND-Zen (^^^Xe), MAJORANA/GERDA (^^Ge), Su- 
perNEMO {^"^Se), SNO-H {^^^Nd), and others [H, with a sensitivity 

Im^^l ~ a few 10-2 eV (14) 
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will probe the IH of neutrino masses. In the case of the normal mass hierarchy |m^^| is much too small in order that 
0:^/3/3-decay will be detected in experiments of the next generation. 

If the Oi//3/3-decay will be observed the measured half-life Ty^'^^^ with experimental error deajp can be converted 

into an "observed effective Majorana mass" |m/3^|°'"' and its error a/jp as 



+ (^)- (15) 



Here, ath is a "theoretical error" of the nuclear matrix element |M''''^|. 

In Fig. [2]we plot the difference of the CP Majorana phases a2i as function of |m^/3|°''^ for the inverted hierarchy of 
neutrino masses and by assuming 0% (blue region), 15% (red region) and 25% uncertainty (orange region) in |to^^|°''*. 
The current experimental errors of neutrino mixing parameters and mass squared differences are taken into account. 
The lightest neutrino mass m-o is assumed to be within a range from to 10 meV and one of the CP violating phases a-i 
(or 0:2) is taken to be arbitrary. We see that if the considered accuracies are achieved it can be possible to determine 
the value of the CP phase difference ai2- However, for cr;3^/|m^^|°''* > 50% it will be difficult, or impossible, to gain 
reliable information about the value of ai2. 

The uncertainty ath in the calculated Oi//3/3 decay NME is a complicated and more involved problem. M'^^ consists 



of the Fermi (F), Gamow-Teller (GT) and tensor (T) parts [29|-l3l| 



^'"=f^) (16) 



Here, g'^/^ {9a) is the quenched (unquenched) axial- vector coupling constant. A/'"" is a function of {g'^^)'^, which 
appears in the Fermi matrix element and also enters in the calculation of the Gamow-Teller and tensor constituents 
due to a consideration of the nucleon weak- magnetism terms [13]. This definition of M'°^ (29l - [3ll | allows to display 
the effects of the uncertainties in g"^-^ and to use the same phase factor G^'' when calculating the Oi//?^-decay rate. 

The treatment of quenching g*^^ is an important source of difference between the calculated Oi//3/3-decay NMEs 
[30I I33I I . Quenching of the axial- vector coupling was introduced to account for the fact that the calculated strengths 
of the GamowTeller /3-decay transitions to individual final states are significantly larger than the experimental ones. 
Formally this is accomplished by replacing the true vacuum value of the coupling constant gA ~ 1.269 by a quenched 
value g"^^ = 1.0. It is not clear whether a similar phenomenon exists for other multipoles besides J = 1^. 

Different nuclear structure methods have been used for the calculation of M'*^"^, in particular, interacting shell model 



(ISM) [3J, |35|, quasiparticle random phase approximation (QRPA) |29l - l3ll l36l]. projected Hartree-Fock Bogoliubov 
approach (PHFB, PQQ2 parametrization) ) [S^, energy density functional method (EDF) 38] and interacting boson 
model (IBM) By assuming an unquenched gA the ISM values of the NMEs are about a factor 2-3 smaller than 
the NMEs of other approaches (see Table 3 of 19]). The results of the QRPA, IBM, EDF and PHB approaches differ 
by a factor less than two. Their results agree rather well with each other in the case of the Oi^/3/3-decay of ^"^"Te. 

A detailed study of uncertainties in the calculated Oz^/3/3-decay NMEs was performed within the QRPA approach 
[29l - l3l1 . IsgI] . The average matrix element (Af'"") (averaged over different NN potentials, choices of the single particle 
space, variants of the QRPA approach) was evaluated, as well as its variance a : 



N 

2 



^5:(A/r-(AO)'. (17) 



" N- 

Further progress was achieved by performing a self-consistent calculation of the NMEs in which the pairing and 
residual interactions as well as the two-nucleon short-range correlations were derived from the same modern nucleon- 
nucleon potentials, namely, from the charge-dependent Bonn potential (CD-Bonn) and the Argonne V18 potential 
[31I ]. The particle-particle strength was adjusted to the 2i//3/?-decay half-life eliminating one of the main reasons for 
variability of the calculated Af '""^ within the QRPA-like methods [2^, . 

Here, we update the calculation of the average {M' ^) and its variance a for the Oz^/3/3-decay of ^^Ge and ^^^Xe. 
The recommended half-life value of 2i^/3/3-decay of ^^Ge [40| and the recently measured half-life of the 2i//3/3-decay of 
^^^Xe [U, H^l were considered. The calculations were performed for the CD-Bonn and Argonne potentials, 3 different 
sizes of the model space, a quenched and unquenched axial vector coupling constant and two variants of the QRPA - 
the standard QRPA and the renormalized QRPA (RQRPA). Following values were found: 

{M'"") = 4.78, a = 0.74 C^Ge) 

= 2.33, a = 0.36 C^'^Xe). (18) 
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FIG. 3. (Color online) The spread of the NMEs for the 0!//3/3-decay of ^^Ge and "''Xe calculated within QRPA and RQRPA 
by assuming Argonne (A) and CD-Bonn (B) NN interactions, 3 different sizes of the model space and quenched {g'^/^ ~ 1.00) 
and unquenched {g''/^ = gA ~ 1.269) weak coupling constant gA- 



The nuclear radius was assumed to be R=1.2 A^^^ fm. We see that the variance a is about 15% of the full NME 
M'"^ Of course, these results are valid only for the QRPA approach and the considered averaging scheme. 

The calculated sets of 24 NMEs for ^*^Ge and of 24 NMEs for ^^exe are presented in Fig. H The results do 
not depend much on the size of the model space and on the type of the NN interaction. For a quenched weak 
coupling constant the NMEs are significantly smaller than those for an unquenched gA mostly due to the factor 
(1.00/1.269) = 0.62 entering in the definition of Af'"'^ in Eq. (|16p. The RQRPA predictions are suppressed when 
compared to the QRPA ones, especially for ""^Ge. 

It goes without saying that further progress in the calculation of the Qv (3 fi-Aec&y NMEs is required. Thanks to 
the theoretical efforts made over the last years the disagreement among the different NMEs is now much less severe 
than it was some years ago. Currently the main issue of interest is that there exists significant disagreement of 
the ISM results with those of other approaches and the problem of quenching the axial-vector coupling constant. 
The uncertainty associated with the calculation of the Oi^/3/3-decay NMEs can be reduced by suitably chosen nuclear 
probes. Complementary experimental information from related processes like the 2^'/3/3-decay [40| . charge- exchange 
[44| and particle transfer reactions [45| is also very important. The differences among results of various nuclear 
structure approaches could be understood by performing an anatomy of the Oi^/3/?-decay NME [s^, [H,!!^ . A recent 
development in the field is encouraging. There is a reason to believe that the uncertainty in the Oj^/3/3-decay will be 
further reduced. 



IV. SUMMARY 



The possible establishment of CP-violation in the lepton sector is one of the most challenging problems of neutrino 
and astrophysics. Studies of Oi^/3/3-decay driven by Majorana neutrinos can lead to insights into CP violation in this 
sector. In view of recent measurement of the smallest neutrino mixing angle ^13 in the Double Chooz, Daya Bay and 
RENO experiments we revisited the problem of determining the Majorana CP phases via the additional observation 
of the Oi^/J/J-decay. 
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Both cases of the normal and inverted hierarchy of neutrino masses were discussed. It was shown that in the case 
of the NH the determination of one of the Majorana CP phases would be possible only by knowledge of both the 
absolute value of the effective Majorana mass |m/3^| and of the lightest neutrino mass mi. This task cannot be solved 
by any of the planned or prepared neutrino experiments, only within some partcle physics models which allow for a 
prediction of neutrino masses. Further, it was also found that for some values of mi the minimal value of |m^^| is 
realized in the case of CP violation. 

The case of the IH of neutrino masses offers different possibilities. Future cosmological measurements have the 
potential to constrain the lightest neutrino mass mo to values below 10 meV. The difference of Majorana phases 0:21 
depends very weakly on mg for these low values and it can be determined by an accurate measurement leading to 
|m/3^|. For this purpose the Oi^/3/3-decay NME needs to be evaluated with an uncertainty less than 30%. This is a 
formidable task, which might be achieved at some point in time due to further developments in the fields of nuclear 
structure and many-body physics also linked to a further increase of computer power. 
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